A numerical investigation of a piezoelectric surface acoustic wave interaction with a 

one-dimensional channel 



O 

o 



CD 



I 

O 

o 



> 
00 

m 

in 
o 



I 

O 

o 



S. Rahman^0 M. Kataoka^, C. H. W. Barnes^, and H. P. Langtangen^ 
^ Cavendish Laboratory, Madingley Road, Cambridge, CBS OHE, United Kingdom and 
'^Simula Research Laboratory, Martin Linges v 17, Fornebu P.O.Box 134, 1325 Lysaker, Norway 

(Dated: February 5, 2008) 

We investigate the propagation of a piezoelectric surface acoustic wave (SAW) across a 
GaAs/AlxGai_x As heterostructure surface, on which there is fixed a metallic split-gate. Our 
method is based on a finite element formulation of the underlying equations of motion, and is 
performed in three-dimensions fully incorporating the geometry and material composition of the 
substrate and gates. We demonstrate attenuation of the SAW amplitude as a result of the presence 
of both mechanical and electrical gates on the surface. We show that the incorporation of a simple 
model for the screening by the two-dimensional electron gas (2DEG), results in a total electric po- 
tential modulation that suggests a mechanism for the capture and release of electrons by the SAW. 
Our simulations suggest the absence of any significant turbulence in the SAW motion which could 
hamper the operation of SAW based quantum devices of a more complex geometry. 

PACS numbers: 85.35.Gv, 73.23.-b, 77.65. Dq 



I. INTRODUCTION 

Surface acoustic waves (SAWs) are widely used in 
microwave circuit components, such as filters and 
resonators. ■'^ In condensed matter physics research, SAWs 
have been a useful tool in probing electronic structure, for 
example, in thin metal films, ^ and quantum Hall liquids,'^ 
for a number of decades. In 1996, Shilton et al. realized 
a device which carries a quantized number of electrons 
through a quasi-one-dimensional channel (QIDC) using 
SAWsji Such devices are currently being developed for 
metrological applications and for quantum logic circuitsi^ 

Despite the enthusiasm for using SAWs in such appli- 
cations, experimental progress in this direction has been 
slow. One of the main reasons for this is that we lack a 
rigorous understanding of the dynamics of the SAW prop- 
agation in complicated device structures, and are there- 
fore unable to precisely design or simulate the devices. 
Earlier analytical work on SAW single-electron transport 
(SET) devicesSiLS required crude approximations, espe- 
cially regarding the effect of surface gates, where a realis- 
tic two-dimensional model of a split-gate device had not 
been achieved. 

In this paper, we present the results of a nu- 
merical investigation of the dynamics at the depth 
of a two-dimensional electron gas (2DEG) in a 
GaAs/AlxGai_xAs heterostructure, resulting from the 
propagation of a piezoelectric SAW through a realistic 
split-gate. The numerical method implemented,^ fully 
takes into account the geometry and material composi- 
tion of the device, as well as the full two-way coupling 
between the electrical and mechanical fields. Our simu- 
lations are performed in the strongly screened, low SAW 
power, low barrier height, regime, as this allows us to 
implement a simple model for the 2DEG and allows as- 
sumptions on the SAW amplitude and split-gate barrier 
height. 

We use our solutions to discuss the effects due to the 



presence of mechanical and electrical surface gates. Al- 
though in this paper we restrict our attention to a split- 
gate device, our method is in principle applicable to a 
device of any geometry. 

This paper proceeds as follows. In Section |nl we de- 
scribe the details of the underlying theory of SAW prop- 
agation in a gated device. We discuss the earlier work 
of Aizin and GumbaSil pointing out the differences be- 
tween our formalisms and solutions. In Section IIIII we 
describe the method behind our numerical strategy. The 
results for an un-gated device and a split-gate device are 
presented in Sections IIV Al and IIV Bl respectively. 



II. SINGLE-ELECTRON SAW DEVICES 

The quantization of acoustoelectric current was first 
demonstrated experimentally by Shilton et almii& Fig- 
ure^a) shows the experimental setup used. It consists of 
a GaAs/Aljf Gai_xAs heterostructure containing a two- 
dimensional electron gas (2DEG) that is formed into a 
mesa by wet etching. A QIDC is formed when a nega- 
tive potential is applied to a split metallic surface gate. A 
large negative voltage applied to the gates induces a nar- 
row depleted constriction between the two 2DEG regions. 
Application of a microwave signal to an interdigitated 
transducer excites electromechanical waves through the 
piezoelectric effect which include both SAWs and bulk 
waves. 

The electric component of the SAW captures electrons 
from the 2DEG in its potential minima transporting 
them through the constriction produced by the split-gate 
potential. The current measured from the drain region 
exhibits a step-like behavior as shown in Fig.^b). The 
values of the current plateaus are quantised with 
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where / is the current, e is the electron charge, n is the 
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FIG. 1: (a) A schematic diagram sliowing the experimental 
setup in acoustic charge transport experiments. A transducer 
on the left excites a SAW wave that propagates towards the 
spit-gate on the right, (b) The acoustoelectric current versus 
split-gate voltage. The current exhibits a step-like behavior 
as the gate voltage is varied. 



number of electrons transported, and / is the SAW fre- 
quency. The case with n — 1 involves the transport of a 
single electron per cycle through the QIDC. 

The original explanation for the quantization of acous- 
toelectric current given by Shilton et al.'^'^'^ asserts that 
the combination of the SAW electrostatic potential and 
the split-gate potential produces a travelling quantum 
dot that transports a fixed number of electrons from one 
side of the constriction to the other. This is plausible 
as the length of the depleted region formed by the split- 
gate (~ 1.5 fiui) is a small fraction greater than SAW 
wavelength (~ 1 /im). Several electrons are believed 
to be initially trapped in each SAW minimum but as it 
passes through the constriction formed by the split-gate, 
its physical dimensions decrease and Coulomb repulsion 
between electrons forces electrons out, thus determining 
the final number of electrons remaining in the dot. This 
explanation, although qualitatively satisfactory, was not 
at the time, supported by a detailed theoretical analysis 
involving an accurate model for the SAW and split-gate, 
which we aim to develop in this paper. 

One approach for such an analysis, is to solve the equa- 
tions of motion in the GaAs/Al^Gai^xAs media, taking 
into account the complex geometry of the split-gate, the 
material composition and external electric fields, to de- 
termine the evolution of the SAW potential through the 
device. The resulting electric potential could be used as 
a part of a quantum-mechanical treatment to describe 
the evolution of the electronic states. 

However, this approach is a non-trivial one, as the 
equations of motion in a piezoelectric material such as 
AlxGai_xAs, consist of a set of four coupled partial 
differential equations, second order in time and space. 
Moreover, as the SAW passes through the split-gate po- 
tential, initially ^ 30 electrons are captured from the 
2DEG, this number reduces as the SAW quantum dot 
is dragged into the increasing electric field in the con- 
striction. This problem then involves a solution to the 
time-dependent many-particle Schrodinger equation in- 
cluding spin. To perform this calculation exactly is be- 
yond the scope or intention of this paper. The focus of 



this paper is the solution of the equations of motion in the 
GaAs/AlxGai-jjfAs medium and, using a simple model 
for the 2DEG and, ignoring any self consistent effects, 
the calculation of the SAW potential as it traverses the 
split-gate. The equations of motion in a heterogeneous 
piezoelectric material are^*ii 

^ dt^ ~ dx,'''^''' dxk dxj^""'' dxk' ^' 

dxi dxj dxi dxj ' 

where Ui is a three-component vector representing the 
displacement of the material at each point, is the elec- 
tric potential, c^^.; are components of the elastic ten- 
sor, (= e"^) are components of the material dielec- 
tric (the superscripts E and S indicate that the values 
were determined under constant electric field and strain, 
respectively), g is the mass density and e^fe are compo- 
nents of the piezoelectric tensor representing the coupling 
strength between the electric and mechanical fields. The 
Eqs. Q can be regarded as a system of wave equations 
with a load term due to the electric field, and Eq. ^ can 
be regarded as a Poisson equation with a source term due 
to the mechanical deformations. 

The equations are solved with traction-free boundary 
conditions. For example cjiz = for the surface normal to 
the z axis, where try are components of the stress tensor, 
defined by 

CTjj = -CktjEk + cfjki£kh (4) 

where E^ = and e., = + The electro- 

static boundary conditions require the prescription of the 
normal component of the electric displacement at the free 
surface of the medium. The electric displacement vector 
is defined by 

D., = efjEj+eijk£3k- (5) 

Aizin et al& followed this approach providing a closed 
form analytic formula for the combined potential of 
the SAW and split-gate. The calculation involved the 
simplification of Eqs. ^ and and of the descrip- 
tion of the physical system, so that analytical solutions 
were tractable. In particular, their approach was two- 
dimensional, taking into account the direction of SAW 
propagation and the direction into the bulk of the het- 
erostructure only. Clearly, this analysis cannot take into 
account the geometry of the surface gates and therefore 
the 'split' nature of the split-gate. In fact, a split-gate 
was approximated to be a single bar gate. Another sim- 
plification was the decoupling of the mechanical motion 
from the electric fields by setting the piezoelectric term in 
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Eq. |5J to zero. This simplification was based on the fact 
that Bijk <C Cijki for AlxGai_xAs. In practice, however, 
externally applied voltages (and hence electric fields) may 
be large enough that they contribute to the mechanical 
strains and hence the observed potential in the vicinity 
of the gates.^ Thus, the term with eijk in Eq. (O cannot 
always be neglected. Their calculations also ignore differ- 
ences in the elastic constants and mass density between 
the gates and the heterostructure substrate. The gates 
are often fabricated from a combination of aluminium 
(Al), gold (Ag) or titanium (Ti). In fact, this approach 
neglects the mechanical presence of the gates at all: the 
SAW potential (j)sAW is effectively assumed to have a co- 
sine form i.e. (f>sAW ^ cos{kx — cut) and is added to the 
split-gate potential (j)sG which is of a quadratic form, to 
give the total electric potential 0. In practice however, 
the mechanical presence of surface gates could cause the 
SAW potential to be damped, reflected or diffracted in 
the vicinity of the gates. The presence of surface gate 
screening has been utilized in SAW-based photolumines- 
cence experiments»i2i^ 

However, using the potential </>, Aizin and Gumbs were 
able to show that if tunnelling of electrons from the quan- 
tum dot to the source 2DEG was allowed, they could ex- 
plain the transition in the acoustoelectric current from 
the point where no charge is transported, to that where 
a single electron is transportedi^ In effect, they provided 
an explanation of the first current plateau. In subsequent 
theoretical work,-^ Gumbs et al. showed, using the same 
potential, that the second quantized plateau in the cur- 
rent as a function of the gate voltage or SAW power can 
be explained by the effect of both the Coulomb blockade 
in the quantum dot and the backward tunnelling into the 
QIDC. 

The analytic approach by Aizin and Gumbs was a con- 
venient one, and similar approaches for the SAW and 
gate potential have been adopted in other theoretical 
worksili^ 

Our recent numerical investigation,^ based on the finite 
element formulation of Eqs. Q and Q demonstrated 
both the effect of the mechanical presence of a gate on the 
SAW electric potential and the presence of a large static 
electric potential on the mechanical strains in the mate- 
rial. The procedure includes the full two-way coupling 
between the electric and mechanical fields. The advan- 
tage of the approach is that we can handle the geometry 
and material composition of the split-gate, or any shape 
of gate pattern, easily. In particular, the SAW wavefront 
was shown to be damped and scattered somewhat, af- 
ter propagating through a single gate. Such effects could 
prevent the proper functioning of SAW devices, operating 
in the quantum regime. However, the intention of that 
paper was to demonstrate that the numerical method is 
capable of reproducing the fully coupled elasto-electric 
dynamics of the physical system. In order to exagger- 
ate the effects of the electrical and mechanical coupling, 
the investigation used fictitious materials for the gate, 
which had substantially different physical properties to 



Al, Ti or Ag, which are frequently used in SAW based 
SET experiments. 

In this paper, we apply our numerical solution strategy 
to a split-gate device with realistic parameters for the 
gate geometry and material composition. We also include 
a simple model for the 2DEG, when a voltage is applied 
to the spht-gate. 



III. METHOD 

Figure |21 illustrates the geometry of the device to be 
modelled. The device consists of a GaAs/AlxGai_xAs 
{X — 0.3) heterostructure with a Ti/Al split-gate placed 
on the surface. The sample has dimensions 12 /im x 
3 /im X 8 /im in the x, y and z directions, respectively, 
with z = zq (zq = 8 /im) at the top surface and the het- 
erostructure occupying the region shown in Fig. |21 From 
the top surface down, there is a 10 nm GaAs capped layer, 
an 80 nm Alj^Gai-jc As layer, and the remainder being 
GaAs. We are interested in calculating the dynamic ef- 
fects of the SAW propagating through the constriction 
and therefore we ignore static charge distributions due 
to doping in the AlxGai_xAs layer. The thickness (or 
height) of the Ti and Al components of the gates are 
chosen to be 20 nm and 40 nm, respectively. The con- 
striction formed by the split-gate has dimensions 0.7 //m 
and 1.0 /im in the x and y directions, respectively. The 
gate is centered along the x axis at x = Xg ~ 9850 nm. 

To obtain an accurate description of the potential land- 
scape in the device (without a SAW) and in particu- 
lar at the 2DEG, a three-dimensional Hartree or density 
functional calculation would be necessary^ incorporat- 
ing suitable boundary conditions for the split-gate and 
charge from donor impurities. Moreover, to incorporate 
the time-varying polarization charge induced from the 
SAW, this calculation would need to be performed at 
every time-step. Such a task would require enormous 
computational resources. Instead, a simple model for 
the split-gate induced potential landscape is implemented 
here, based on our experimental observations. In the ex- 
periments, the split-gate potentials are applied relative 
to the 2DEG. The 2DEG is believed to be depleted a few 
hundred nano-meters laterally away from the split-gate. 
Hence, in this model, when a gate voltage is applied, 
the 2DEG is assumed to be depleted 300 nm laterally 
from the split-gate, and is modelled as a metal sheet, 
with a potential difference applied between the split-gate 
and the metal sheet (with the potential at the 2DEG set 
to V), in the numerical solution of the Eq. For 
the electrical boundary conditions at the free surface, 
we choose a Neumann-type condition on the un-gated 
regions, implemented by setting = in order that 
overall charge neutrality in the device is satisfied>i2iiS 
Dirichlet-type conditions are unsuitable for the un-gated 
regions as they would compromise the SAW potential. 
The mechanical boundary conditions implemented are 
the 'traction-free' conditions defined earlier. 
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4.416 


n/a 
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2.626 


n/a 
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TABLE I: Physical properties of AlxGai-xAs^ A^ and T^. The permittivity of vacuum is eo. 



Time-dependent boundary 
condtion applied here 



2DEG 



(0, 0, 0) 









T 

(X 


, 0, z„) 



FIG. 2: A schematic diagram showing the device geometry. A 
split-gate, composed of a Ti/Al alloy, is placed on the surface 
(2 = zO) of the GaAs/AlxGai_xAs substrate. The SAW, 
generated by the three membranes separated 1.0 fim apart, 
propagates along the x-axis. The xo, yo and zo parameters are 
chosen to be 12000 nm, 3000 nm and 8000 nm respectively. 
The depletion parameter d is 300 nm. The gate is centered 
along the x axis at x = Xg = 9850 nm. 



In order to excite the SAWs, we apply a time- 
dependent boundary condition of the form 



Uz = Asin(27r/<), x ~ const, z — zq, (6) 

where / is the frequency and A is the amplitude. Here, 
/ is chosen to be 2.7 GHz in order to satisfy the re- 
lation vsAW = /A, where vsaw is the SAW velocity 
(~ 2.7 X 10^ ms~^) and A is the period of the transducer 
fingers. The amplitude A is chosen such that the SAW 
has the amplitude of ^ 20 mV at a depth of 100 nm, 
and corresponds to the low SAW power regime of the ex- 
periments. The oscillation is applied at three values of x 
separated by 1.0 /im as shown in Fig. [3 to increase the 
SAW amplitude in the x direction relative to z direction: 
the SAW wave-fronts generated by each membrane add 
more constructively along the x-axis than the z-axis. 

For AlxGai_xAs, the only non-zero components of 
the piezoelectric tensor have the value 614. Also, the 
non- vanishing components of the permittivity tensor are 
£11 = £22 — £33 — Cs, and the non- vanishing components 
of the elastic tensor (not written out above) are Cxxxx — 
^yyyy — ^zzzz — Cn, Cxxyy — (^yyzz — ^zzxx — C12, and 
Cxyxy = Cyzyz = Czxzx = C44 . All Other non-zcro compo- 
nents of the elastic tensor can be determined by applying 
its symmetry properties ci^u = Cjiki = Cyifc = cuij. The 
values of these constants for the relevant materials are 
given in Table I. The dielectric constant and ei4 are not 



required for the Ti and Al as they are subject to Dirich- 
let boundary conditions on all external surfaces. In our 
simulations, we follow the convention in SAW based SET 
experiments where the crystal orientation is such that the 
x direction is aligned to [Oil] direction, and z direction 
to [100] direction.20 

The finite element method was chosen for its 
proven ability in handling geometrically complicated 
domains The basic idea of the finite clement method 
is to approximate the unknown fields, for example (j) in 
the Poisson equation above, by a linear combination of 
basis functions Ni, cj) k cf) — J2^j^i^j4'j: then insert (j) 
in the Poisson equation, and demand the residual to be 
orthogonal to the space spanned by {TVi, . . . , A^„}. We 
utilize finite elements in the spatial part of the problem 
and use a second order finite difference discretisation in 
time for computational speed so Eq. (|2Jl is approximated 
by 



,,e-i 



2u' 



At2 



RHS', 



(7) 



where RHS is a finite element approximation of the elas- 
tic stress term and the mechanical field load. The su- 
perscript £ represents the time level. We use eight-noded 
brick elements corresponding to linear basis functions Ni, 
resulting in an overall spatial and temporal convergence 
rate for the error, of 2.0. An operator splitting strategy 
is employed to split the coupled Ui-ip problemi^ Inheri- 
tance and polymorphism principles from object-oriented 
programming, 22 are used to incorporate software com- 
ponents from the Diffpack library, for solving the Pois- 
son and the elasticity equations, respectively, thus max- 
imizing efficiency in the programming and verification. 
The numerical formulation and verification is described 
in more detail in Ref. |2 



IV. RESULTS AND DISCUSSION 

A. SAW through an un-gated surface 

We first demonstrate that our method for exciting 
acoustoelectric vibrations using three vertically oscillat- 
ing membranes, excites the SAW modes with the re- 
quired wavelengths and velocity, by performing a simu- 
lation without any surface gates and without the 2DEG. 
Therefore, it is not necessary to set the total electric po- 
tential (/) to V at the 2DEG in this simulation. (Also, 
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the absence of depletion would result in the SAW being 
invisible everywhere at that depth.) Figure |31 shows the 
time development of the electric potential as a function of 
distance along the propagation direction (+x direction). 
The curves were extracted at a depth of 100 nm (where 
the interesting physics of acoustoelectric charge trans- 
port take place). The time-dependent boundary condi- 
tion given by Eq. ® is applied on the left-hand-side of 
Fig. 13 with the closest membrane at cc = 6600 nm (not 
shown). We see that two transient peaks, clearly distin- 
guishable from all the other peaks, are initially formed, 
and are followed by a more consistent set of peaks which 
also propagate from left to right of the plots. The tran- 
sient peaks, which are a common feature of these nu- 
merical simulations, are a consequence of the transition 
from a flat lattice i.e. zero displacement everywhere at 
(i — 0) to one with SAWs {t > 0). Moreover, as three 
membranes are used to generate the SAW, it takes three 
periods of the SAW to establish a consistent set of peaks. 
In experiments the transients occur too, but many thou- 
sands of SAW wavefronts pass through the split-gate, so 
the initial transient peaks are actually insignificant. In 
the following numerical experiments, the transients are 
included in the analysis as they provide some insight of 
differing SAW amplitudes, in a single simulation run. 

From Fig. |31 we see that in this particular run, the 
largest amplitude is ~ 20 mV and the SAW resem- 
bles the plane wave shape with the expected velocity 
2770 ± 20 ms~^, and a wavelength of 1.0 /im. In our 
previous numerical investigations^^ we showed that the 
solution method exhibited the characteristic decay into 
the bulk as well as a phase difference between the lat- 
eral and vertical displacements Ux and Mz, respectively!^ 
From Fig. O we see that despite the discontinuities in 
the material parameters due to the the presence of the 
heterostructure, we are able to produce coherent waves. 
This is consistent with the fact that the electrical and me- 
chanical properties of Al^Gai-xAs as shown in Table I, 
differ by less than 1 %. 



B. SAW through a split-gate 

The mechanical effects of the gates are investigated 
first, by performing a simulation of a SAW through 
a split-gate and without applying any gate voltages. 
Dirichlet boundary conditions are therefore not applied 
to the gates in this case. Moreover, since here we are in- 
terested only in the effect of the purely mechanical pres- 
ence of the gates, we can assume that they have dielec- 
tric properties identical to that of the substrate below 
it. This allows us to avoid divergences in the solution of 
Eq. ||2Jl due to absence of the dielectric parameter , for 
the metals. Figures0Ja) and^fb) compare the SAW po- 
tential 4> with and without the mechanical surface gates, 
directly below the surface gate and at the center of the 
constriction, respectively. Directly below the gate, the 
SAW amplitude is seen to be reduced by ~ 30 %, al- 
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FIG. 3: Time development of SAW induced electric potential 
(/!> as a function of the a;-axis (without gates on the surface and 
without the 2DEG). For clarity, the inset shows the amplitude 
of the SAW is ~ 20 mV. 



though the amplitude recovers up to at least 90 % of 
its original value after traversing the gate. At the center 
of the constriction, the difference of amplitudes between 
the gated and un-gated devices is less than one per-cent 
in these simulations. The results of Fig. ^ are consistent 
with the fact that electromechanical energy carried by 
the SAW is much greater than that could be stored by 
the gates; the size of the SAW into the depth is over one 
micron while the height of the gate is typically ~ 0.05 /im. 

Having established the mechanical effects of the Ti/Al 
split-gates, we now apply an electric potential to the sur- 
face gates such that the barrier height due to the split- 
gate, after ignoring static charges on the surface and 
donor levels, is ~ 20 mV at 100 nm below the surface, 
corresponding to the low barrier height and short con- 
striction, regime in experiments. From the initial time- 
level of the simulation at t = to the final at t = 2.5 ns, 
several SAW wavelengths pass through the depleted re- 
gion. We observe a pattern from the total SAW and gate 
potential (j), resembling that of an 'electron pump'^SLSi 
which is repeated with the period of the SAW. Figure 
actually includes three cycles of this modulation. For 
the third modulation starting at < = 2.24 ns, a local 
minimum extends between the source 2DEG and the de- 
pleted region. We would expect electrons from the source 
2DEG to relax into the minimum. As t increases, a SAW 
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FIG. 4: The SAW potential <p mt t = 2.0 ns, (a) directly under 
the gates (y = 500 nm) (b) in between the gates at the center 
of the split-gate constriction (y — 1500 nm). The curves 
were taken 100 nm below the surface. The SAW amplitude is 
reduced directly below the gate as much as ~ 30 % but less 
than 1 % at the center of the constriction. 



maximum enters the depleted region forming a potential 
barrier against electrons escaping or tunnelling backward 
into the source 2DEG while the SAW maximum in front 
acts against escaping or tunnelling forward into the drain 
2DEG. The electrons are thus confined in the SAW min- 
imum and are transported along with it, over the poten- 
tial barrier formed by the split-gate, until the forward 
SAW maximum leaves the depleted region and becomes 
screened by the drain 2DEG. 

Figure [S] suggests that the first two pump motions, 
which are caused by the transients in the SAW and thus 
have a smaller amplitude, have a lower probability to 
transport electrons through the channel, as the potential 
minima that extend from the source 2DEG to the de- 
pleted region are smaller, producing overall smaller min- 
ima when combined with the gate potential. Therefore, 
fewer states will be available in the dot for the electrons to 
occupy. For these transients, the rear tunnel barrier are 
also smaller and therefore electrons have a higher prob- 
ability to escape to the source 2DEG. This is consistent 
with Gumbs' investigation^i which showed the acousto- 
electric current increase as the ratio of the SAW potential 
amplitude to the height of the gate induced barrier, is in- 
creased, despite a different reasoning based on analytic 
models for the SAW and split-gate potentials. 

Experiments have been performed in low SAW power 
and low barrier height regime, which utilized SAWs in 
conjunction which a quantum dot to effect a pumping 
modulation of the total electric potential. 

Figure |B1 shows the curve plots of the total electric po- 
tential <j) for times t between 2.24 ns and 2.34 ns when 
a SAW maximum enters the depleted region. It is clear 
that the rear potential barrier becomes wider and taller 
with time, thus reducing the probability of electrons es- 
caping or tunnelling backwards into the 2DEG. These 
results suggest that the number of electrons transported 
through the dot is determined early on in the capture 
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FIG. 5: Time development of SAW induced electric potential 
as a function of the x axis. The simulation includes the 
screening by the 2DEG, in the source and drain regions. 
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FIG. 6: The time development of a rear potential barrier 
showing the mechanism by which electrons from the source 
2DEG could be captured. 



process. 

In much of the previous theoretical workjSiS the in- 
crease in width and height of the rear potential bar- 
rier has not formed a significant part of the analysis. 
However, Flensberg et alm^ demonstrated that the rapid 
change of the SAW potential barrier at the entrance 
to the channel results in a rapid reduction of tunnel- 
coupling between the source 2DEG and the SAW min- 
imum inducing non-adiabaticity in the travelling dot. 
Flensberg et al. then showed that the non-adiabaticity 
sets a limit on the accuracy of the quantization plateaus. 

In the model of Robinson and Barnesi^ the tunnel bar- 
rier decreases with time, and the number of electrons in 
the dot is determined when the SAW minima reach the 
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(c)f = 2.375 ns (d) ? = 2.500 ns 



FIG. 7; Surface plots showing the time development of the 
SAW induced electric potential <j) in the x-y plane. 

point of maximum gradient of the split-gate potential 
barrier. However, their work was based in the high SAW 
power and high barrier height, regime and cosine and 
gaussian models for the SAW and split-gate, respectively, 
were implemented. 

Figure 13 shows surface plots of the electric potential 
on a two-dimensional surface on the x-y plane, 100 nm 
below the surface at a sequence of times t providing a 
more vivid illustration of the pumping mechanism. In 
particular, aX t = 2.35 ns a quantum well begins to form 
on the left hand side of the split-gate as a SAW mini- 
mum enters the constriction. At t = 2.375 ns the SAW 
minimum is located in between the gates, forming a well 
defined quantum well of a circular geometry (although 
the geometry will vary from circular to elliptical depend- 
ing on the split-gate potential), ki t — 2.5 ns, the SAW 
minimum exits the constriction, and the well begins to 
disappear. 

FigureslH^a) andlH^b) show curve-plots of the total elec- 
tric potential parallel to the y-axis, through the center 
of the constriction. From Fig.lHfa), one can see the oscil- 
lations of the total electric potential as the SAW prop- 
agates through it. The amplitude of the oscillations at 
the center of the constriction is ^ 20 mV, about the same 
amplitude as that of the bare SAW (i.e., without gates), 
whereas the SAW amplitude below the gate is ~ 5 mV 
and therefore has been reduced by ~ 75 %. This is due 
to the screening of the SAW potential by the gates, the 
presence of mechanical gates as discussed earlier, and to a 
much lesser extent, due to the mechanical strains caused 
by the gate voltage affecting the SAW motioni^i^ These 
figures also show the a change in the SAW frequency 
under the gates, in the form of 'splitting' of the peaks, 
where labelled. In particular, the SAW oscillates at a 
higher frequency under gates. 

The peak splitting or increase in the frequency of the 
SAW under the gates is likely to be caused by the dif- 
ference in elasticity constants of the gates and substrate. 
In particular, the elasticity constants of the gates are 
greater that of the substrate. The gates are therefore 
much 'stiffer' than the substrate, and have a higher natu- 
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FIG. 8: (a) Time development of the total electric potential 
as a function of the y coordinate, (b) The electric potential 
cj> directly under the gate and in between the gate. 



ral frequency of vibration. The vibration of surface gates 
impart additional vibrational components to the under- 
lying SAW. A more detailed understanding of the mech- 
anism behind the peak splitting would involve investiga- 
tion of the lateral and vertical components of the dis- 
placements under a gate (and perhaps experiments with 
different gate materials), and therefore a major digres- 
sion from the scope and intention of this paper but a 
possible topic of future work. 

The simulated changes in amplitude and frequency of 
SAW would be undesirable for electron transport exper- 
iments as they may allow both tunnel and Coulomb in- 
teractions between electrons in consecutive SAW minima. 
However, these phenomena seem to be more significant 
for the first three peaks passing through the gate, than 
for the last peak. The last peak has a larger amplitude 
and thus greater energy, and is therefore more resistant 
to these effects. 

Finally, from Fig. |H1 one can see the absence of 
any significant turbulence in the SAW motion along 
the y axis, which may have affected systems of one- 
dimensional channels in parallel which are currently un- 
der development , ,30 'pjjjg g^jgp g^jiggests that diffraction 
of the SAW on propagation though our split-gate is be- 
yond the accuracy of our numerical work i.e. extremely 
small. Therefore the only waves in our system are plane 
waves which travel parallel to the lateral edges of the 
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computational domain and are not reflected of them. 

V. CONCLUSION 

We have performed numerical simulations in order to 
investigate the dynamics resulting from the propagation 
of a piezoelectric SAW though a QIDC defined by metal- 
lic split-gates on the surface of a GaAs/Aljf Gai_xAs 
heterostructure. Our simulations were performed in 
three spatial dimensions fully incorporating the mechan- 
ical and electrical parameters of the materials, and were 
based on the strongly screened, low SAW power and short 
constriction regime. We ignored the presence of static 
charges, as we are interested in the dynamical properties 
of the SAW as it propagates through a realistic split-gate 
device. 

We have demonstrated significant 'amplitude reduc- 
tion' of the SAW potential, up to 30 % due to the me- 
chanical presence of the surface gates, as well as up to 
75 % screening of the SAW potential, by the biased split- 
gate device. In addition, the simulations demonstrated 
the recovery of the SAW, after amplitude reduction due 
to the presence of mechanical gates, and also the absence 
of significant damping or screening effects at the center of 



the constriction formed by the split-gate. These effects 
have been demonstrated theoretically for the first time 
for a realistic device, and would be difficult, if possible 
at all, to achieve without a numerical procedure such as 
ours. The results suggest that the coherent propagation 
of SAWs through systems of QlDCs in series or in parallel 
could also be achieved, although for definiteness, further 
simulations may be required. 

Through the incorporation of a simple model for the 
screening by the 2DEG, we demonstrated a total electric 
potential modulation that resembled an electron pump 
and provided a simple model for the capture of electrons 
by the SAW, from the source 2DEG, and the release of 
the electrons to the drain 2DEG. 
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